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Abstract. - We investigate the transport of a single excitation through a chain of weakly coupled 
subunits. At both ends the chain is exposed to baths which are incorporated by means of a master 
equation in Lindblad form. This master equation is solved by the use of stochastic unraveling in 
order to obtain excitation profile and current in the steady state. Completely diffusive transport is 
found for a range of model parameters, whereas signatures of ballistic behavior are observed outside 
this range. In the diffusive regime the conductivity is rather independent from the strength of the 
bath coupling and quantitatively agrees with the diffusion coefficient which has been derived from 
an investigation of the same model without baths. Also the ballistic behavior in the non-diffusive 
regime is in accord with results from this alternative approach. 



> 

oo 

o 

On 
O 



X 



There essentially exist two direct approaches to the in- 
vestigation of gradient-driven transport phenomena such 
as, e.g., heat conduction: (i.) A closed scenario is consid- 
ered, where transport is driven by an internal gradient, 
i|.e., transport is somehow analyzed for the relaxation of 
a spatially non-uniform energy distribution (This relax- 
ation is in accord with a diffusion equation in the case 
of normal transport.) [1-5]; (ii.) An open scenario is con- 
sidered, where transport is induced by an external gra- 
dient, i.e., baths with different temperatures are locally 
coupled to both ends of a system such that a stationary 
non-equilibrium state results (with a finite current and 
a spatially linear energy profile in the normal transport 
case) [3-12]. 

Furthermore, there are two indirect approaches to 
gradient-driven transport: (iii-) The Green-Kubo ap- 
proach [13] which essentially gives a diffusion coefRcicnt 
as an integral over the current auto-correlation function 
corresponding to the spatially uniform equilibrium state; 
(iv.) The Einstein-Herzfeld approach [14] which relates 
the diffusion coefficient to the mean square displacement 
in the equilibrium state. (For a comprehensive fundamen- 
tal review of (iii.) and (iv.) in the context gradient-driven 
transport see Ref. [15] and references therein.) 

However, even though all those approaches are com- 
monly expected to yield equivalent results on transport, 
direct, say, qualitative comparisons for concrete quantum 



models appear to be rare in the literature [5-7, 11]. Quan- 
titative ones are even more rare. For the model addressed 
below the equivalence of (i.) and (iii.) has been investi- 
gated and partially confirmed (in the diffusive regime) in 
Ref. [1]. Thus, in this Letter we focus exclusively on the 
equivalence of (i.) and (ii.), leaving the comparison of, say, 
(i.) and (iv.) as a relevant subject of future work in that 
direction. 

The Letter at hand provides such a comparison for a 
quantum model which describes the transport of a sin- 
gle excitation through a finite chain of weakly interact- 
ing subunits. This chain has already been treated as a 
closed system and is known to exhibit both purely diffu- 
sive and completely ballistic dynamics, depending on the 
model parameters [1,2]. At both ends the chain is now 
exposed to baths which are incorporated by means of a 
master equation in Lindblad form [16, 17] such that the 
efficient numerical method of stochastic unraveling is ap- 
plicable [18,19]. For this open system we also find the 
above regimes of diffusive and ballistic transport. 

In the normal transport case we particularly show the 
following: (i.) The steady state of the open system has a 
constant current and a linear excitation profile; (ii.) The 
resulting conductivity is almost independent from the 
strength of the bath coupling; (iii.) This conductivity 
coincides with the theoretical diffusion coefRcient of the 
closed system, too. The latter two points deserve closer 
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attention, since it generally is difficult to ensure that the 
extracted conductivity is a pure bulk property, especially 
for a finite system [3,4,7]. 

In the ballistic transport case we finally demonstrate 
the sensitive dependence of the conductivity on the bath 
coupling strength. In particular we illustrate that huge, 
say, infinite conductivities can be observed solely in the 
limit of weak bath couplings, e.g., where the resistance 
due to the bath contact is large. 
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Fig. 1: Sketch of the considered open system: a chain with a 
bath at both ends. The chain consists of A'' identical subunits 
which feature a non-degenerate ground state, a wide energy gap 
AE, and a narrow energy band 5e with n equidistant levels. A 
state from the investigated one-excitation subspace is indicated 
(circles). 

Concretely, we investigate an open quantum system ac- 
cording to Fig. 1: It is an one-dimensional structure which 
is connected to baths at both ends. Before the realization 
of these baths is discussed below, we specify the chain it- 
self, i.e., we introduce the closed quantum system at first. 
It consists of TV identical subunits which feature a non- 
degenerate ground state eo = 0, a wide energy gap AE, 
and a narrow energy band 6e with n equidistant levels 
a = AE + i6e/{n - 1). 

In the following the consideration will be focused on 
the invariant zero- and one-excitation subspace which is 
spanned by the basis {|0), |/.t, z)}. In the single state |0) all 
subunits are in their ground state. In the N n states i) 
only the fith subunit is excited to the ith level of its band, 
while all other subunits arc still in their ground state. By 
the use of this notation the local Hamiltonian of the ^th 
subimit may be written as h^^ ~ 127=1 The 
next-neighbor interaction between two adjacent subunits 
jj, and /X -|- 1 is supposed to be 



(1) 



with the overall coupling strength A. The /x-independent 
coefficients Cij are complex, independent, and random 
numbers: their real and imaginary parts are both cho- 
sen corresponding to a Gaussian distribution with mean 
and variance 1/2. Note that only a single realization of 
the Cij (and not some ensemble average over different re- 
alizations) is considered throughout this Letter. The total 
Hamiltonian may be written as H = Hq + V , where Hq is 
the sum of the local Hamiltonians /i^ and V is the sum of 
the next-neighbor interactions Vf^, respectively. 



Of particular interest are the local probabilities (t) for 
finding an excitation of the /ith subunit somewhere to its 
band. These quantities are conveniently expressed as the 
expectation values p^(t) = Tr{p(<)j5^t} of respective oper- 
ators = 1/^' *) (Mi ^1 J where p{t) denotes the density 
matrix of the whole chain. The associated local currents 
j^{t) may be defined as the expectation values of the op- 
erators = I [{p^ - p^+i)/2, fi^], cf. Ref. [1]. 

Even though this system is not meant to represent a 
concrete physical situation, it may be illustrated as a sim- 
plified model for a chain of, say, coupled atoms, molecules, 
or quantum dots. In this case the hopping of the excita- 
tion from one subunit to another corresponds to energy 
transport, especially if A ^ AE. It may also be viewed 
as a model for non-interacting particles on a lattice with 
many orbitals per site. The hopping of the excitation cor- 
responds to transport of particles in this case. A more 
detailed discussion of physical realizations can be found 
in Refs. [1,2]. 

However, for the closed system it has reliably been 
shown that the dynamical behavior of the p^{t) is well 
described by a diffusion equation of the form p^(t) = 
D [p^t„i(i) — 2pi^i{t) + p^+i{t)], if only the two conditions 
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are fulfilled, see Ref. [2]. The pertinent diffusion constant 
is given hy D = 2Trn /Se. Thus, the relaxation time 
of the excitation profile's Fourier component with wave 
length s, i.e., the characteristic time for transport on this 
length scale, is r^j = 1/[2£'(1 — cos27r/s) ]. 

The first condition (2) guarantees that the dynamics on 
the shortest length scale (s = 1) is much slower than a 
generic correlation time as generated by the local parts of 
the Hamiltonian: rc ~ 1/Se. 

The second condition (2) ensures that the dynamics on 
the longest length scale (s = N) is much faster than the 
period of any correlation function as generated by the lo- 
cal parts of the Hamiltonian: T = 2TTTi/5e. When this 
criterion is broken, e.g., by A — > or TV ^ oo, a transition 
towards ballistic transport is caused in the large s-limit. 
In that case the dynamics is still found to be governed by 
the above diffusion-type equation but now with a time- 
dependent rate, i.e., D ^2Dt/T, cf. Ref. [2]. 

Our aim is to compare these findings with the results 
which are obtained in the following treatment of the same 
model as an open system which is exposed to baths. In 
order to incorporate these baths we now postulate that a 
quantum master equation (QME) in Lindblad form holds, 
namely. 



(3) 



The first part on the r.h.s. of this equation describes the 
coherent time evolution of p{t) w.r.t. H . The second part 
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is an incoherent damping term and given by 

T^Pit) = (Apit)Al ^[p{t),AlA,l\ (4) 

i ^ ^ 

with non- negative rates a^, Lindblad operators A^, and 
the anti-commutator [...,.. . Eqs. (3), (4) are the most 
general form of a linear and time-local QME which defines 
a trace- and hermicity-perserving, completely positive dy- 
namical map. This particularly means that any density is 
mapped to a density matrix [16, 17]. 

It is well-known that a strict derivation of a physically 
reasonable Lindblad form from a microscopic bath model 
is somewhat subtle for conduction scenarios [8]. However, 
since numerics eventually turns out to be rather involved 
in that case, here we primarily choose the computation- 
ally least costly bath implementation (within the estab- 
lished framework of open quantum systems) which yields 
a physically acceptable result. Lindblad forms are com- 
putationally preferable, since those allow for the efficient 
numerical method of stochastic unraveling [18,19]. But 
for our below choice of Lindblad operators the underly- 
ing microscopic formulation can in principle be found by 
following the ideas in Ref. [8]. 

Thus, we choose only two left-bath operators 

n 

with rates oo = ai = l'^ but In right-bath operators 

M^ = R+ = 1^-1, i){0\ , ia^+i = R; = , (6) 

i = l,...,n with rates a2i = r+, a2i+i = i+ and 

L~ describe transitions between the ground state |0) and 
a state where the excitation is equally distributed over all 
band levels of the leftmost subunit. Rf and R~ are re- 
sponsible for transitions between the ground state and the 
band state |iV — 1, i) of the rightmost subunit. (Note that 
this contact modelling goes beyond the secular approxima- 
tion and is not in principle conflict with some microscopic 
picture [8]). We further choose r+ = 0, i.e., we suppose 
that the right bath is at zero temperature. The tempera- 
ture of the left bath is finite, cf. Fig. 1. 

The method of stochastic unraveling relies on the 
fact that any Lindblad QME for a density matrix p{t) 
can equivalently be formulated in terms of a stochastic 
Schrodinger equation (SSE) for a wave function \ip{t)), 
see Refs. [18,19]. This fact advantageously allows to deal 
with the same problem in the lower dimensional Hilbert 
space. In order to simulate the process which is defined 
by the SSE we apply the procedure in, e.g., Ref. [7]: 
It consists of deterministic evolutions (w.r.t. an effective 
Hamiltonian HcS = H — i/2 a,; A[ ) and stochastic 
jumps (corresponding to one of the Lindblad operators 
Ai). Each application of this procedure leads to a single 
realization of a so-called "trajectory" \ip{t)). The expec- 
tation value of some observable M in the steady state 



p can be evaluated from such a trajectory by the time 
average Ti{pM} « 1/t /J dti{ip{ti) \M\ ip{ti)) for a sufR- 
cicntly large time interval t [19]. We are primarily inter- 
ested in M ~ and M ~ j^, of course. 

Due to our choice of Lindblad operators (and r+ = 0), 
the procedure drastically simplifies, because exclusively 
jumps into the ground state |0) or into the state L\0) 
can occur, cf. Eqs. (5) and (6). Furthermore, since |0) is 
an eigenstate of Hcff, its deterministic evolutions become 
trivial and feature probabilities Pfj,{t) = and currents 
jp,{t) = 0. Consequently, to the above time average only 
those parts of \(p{t)) contribute, where the deterministic 
propagations of L\0) take place, i.e., the whole problem 
reduces to the evaluation of |^/'(t)) = cxp{—iHcsT)L\0) 
solely. Note that propagations with the length of time 
r appear in \(p{t)) with a frequency of occurrence which 
is given by W{t) = -d/dr (V'(r) | ?/>(t)), cf. Ref. [7]. It is 
therefore possible to derive an exact analytical formula for 
the infinite time average, namely. 
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Tt{pAI} = 








with the constant C = 1/1+ + dr M^(t) t. Note that 
the dependence on the bath rate 1+ appears as the scaling 
factor C solely. Thus, the remaining crucial parameters 
are l~ and r~ . For simplicity, however, we consider the 
special case of a negligibly small l~ (within Hcs), i.e., the 
case of a sufficiently weak coupling of left bath and system. 
(In that case also 1+ becomes small but only affects the 
above constant C. No restriction is made for the ratio 
/l+, e.g., the temperature of the left bath is still finite.) 
Only a single bath parameter remains: = r. 




Fig. 2: Local probabilities (circles) and local currents 
(squares) for (a) A = 0.001, r = 0.0025 and (b) A = 8 • 10"^ 
r = 4 ■ 10"''. Remaining model parameters: N = 5, n — 500, 
and (5e = 0.5 (D = 6280 A^). The resulting conductivities are 
(a) K,/D ^ 0.98 and (b) /t/D ^ 4. 

In practice the above integral (7) is approximated by a 
sum over discrete time steps Ar^ and with a finite upper 
limit Tniax- The use of not too large, i.e., numerically still 
accessible T^ax is possible, because the function W{t) usu- 
ally decreases rapidly with r, e.g., it decays exponentially 
fast. In OUT numerical simulations we use At^ D < 0.01 
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Fig. 3: Conductivity k as a function of the bath rate r for 
A = 0.001 (circles), A = 16 • lO"'' (triangles), A = 10 ■ lO"'^ 
(rhombuses), and A = 8 • 10^'' (squares). Remaining model 
parameters: iV = 5, ?i = 500, and Se = 0.5 (D = 6280 A^). 



and Tmax D > 1000. This choice typically is minimum re- 
quired for a good convergence in terms of which are 
independent from fj,, i.e., with a deviation from their aver- 
age in the order of less than 1%, at least for the following 
Figs. 2 and 3. 

Fig. 2a shows first numerical results for a system with 
N ~ 5, where the remaining model parameters are set 
to well satisfy the both conditions (2) for purely diffusive 
transport in the closed scenario. In the open scenario at 
hand, as expected for normal behavior, we indeed observe 
a strictly linear excitation profile over the whole chain 
without significant effects at the boundaries n = and 

— 1. Moreover, the conductivity k = jfj./{Pfi — Pfi+i) 
turns out to agree almost perfectly with the theoretical 
diffusion coefficient D from Ref. [2], k/D « 0.98. This is 
one of the main results of this work. It already indicates 
that K may be viewed as a bulk property and not as some 
artifact of the bath coupling. So far, k may be extracted 
at any position fi. But later on we will extract k always 
from the middle of the chain. 

However, a detailed investigation of ^(r) for the same set 
of model parameters is presented by the circles in Fig. 3. 
Remarkably, k is found to be independent from r over al- 
most three orders of magnitude. The decrease of k in the 
large r-limit may be understood as an effect which is due 
to the onset of the breakdown of the weak coupling ap- 
proximation that underlies the Lindblad form in general, 
since it starts at those rates which are close to the inverse 
correlation time I/tq ~ Se of the closed system. (Note 
that the underlying excitation profile is still linear ex- 
cept for an abrupt decline to zero at = N — 1.) 

The decrease of n in the small r-limit appears to be 
contra-intuitive but nevertheless is numerically observed 
here. This decrease begins at those rates which are close 
to the inverse relaxation time 1 /ta ^ D cx of the closed 
system. Therefore one may carefully conjecture that the 
onset of the decline is shifted to even smaller r, if only 1 /t^ 



becomes smaller, e.g., if A is decreased. This conjecture is 
confirmed by the triangles in Fig. 3. 

But the further reduction of A eventually leads to the 
violation of the second condition (2), i.e., in the closed sys- 
tem diffusive transport is finally expected to break down 
towards ballistic behavior. As displayed in Fig. 2b, such a 
behavior is found for the open system at hand, too. The 
probability profile is highly non-linear and the conduc- 
tivity K, now extracted in the middle of the chain, is also 
much larger than the theoretical diffusion coefficient D, 
k/D K, A. (Note that the integration of (7) becomes very 
costly, since the function W{t) exhibits a slowly decaying 
long-time tail.) 

Contrary to the above diffusive case, k turns out to be 
extremely r-dependent in the ballistic case, see the squares 
in Fig. 3. (Note that these squares are evaluated for the 
maximum r-interval which is still accessible to our nu- 
merics). K{r) features a distinct peak below the inverse 
correlation period 1/T = 5e/{2irn) of the closed system. 
The decrease of n on the l.h.s. of the peak may be un- 
derstood as the decline to zero which has already been 
observed for r ^ 1/tr. On the r.h.s. of the peak the de- 
crease of K, reaches rather fast a value which is close to D 
and eventually seems to reach this value but to stagnate 
at it. Surprisingly, the ballistic nature of transport does 
not have a significant impact on k here. 

This observation, and all other findings in Fig. 3, allow 
for the following interpretation: Assume that the external 
bath resolves the internal dynamics on a corresponding 
time scale t = 1/r. Then a deviation from the transport 
behavior of the isolated system is to be expected, if the 
bath-induced dynamics is either faster than internal corre- 
lation times (r > I/tc) or slower than internal relaxation 
times (r < I/tb). (This deviation results in a decrease of 
K.) If the bath-induced dynamics proceeds in the possibly 
huge regime in between those time scales, diffusive behav- 
ior becomes visible. However, a ballistic closed system 
typically features a time scale on which the non-decaying 
character of the correlation functions becomes crucial but 
relaxation has not been completed yet. As already ex- 
plained, for our model such a time scale exists, M tr <T 
(i.e., if the second condition (2) is violated), precisely in 
between those two times. Now, if the bath-induced dy- 
namics is tuned by r to the latter time scale, signatures 
of ballistic transport, i.e., a sharply increasing k, appear. 
The larger the latter time regime becomes, the more pro- 
nounced is the maximum of k and its height may probably 
rise without bound according to the degree to which the 
second condition (2) is violated, see Fig. 3. This param- 
eter regime is, however, numerically very costly and can 
with our resources only be explored to the extend dis- 
played in Fig. 3. It is numerically even more costly to 
drive the violation of the second condition (2) by increas- 
ing the length TV rather than decreasing A, which prevents 
us from exploring this regime directly. Nevertheless, based 
on the above interpretation, we expect that the length 
scale dependent transition to ballistic behavior discussed 
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in Ref. [2] should become visible in the present bath sce- 
nario, too. 

To conclude, wc investigated a transport scenario based 
on bath coupling for a model with quantitatively well- 
known transport properties. We found that a bath sce- 
nario may produce the correct bulk conductivity but this 
requires a careful implementation of the bath contacts. 
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